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We present the corrections to the fermion propagator, to second order in the lattice spacing a, in 
1-loop perturbation theory. The fermions are described by the clover action and for the gluons 
we use a 3-parameter family of Symanzik improved actions. Our calculation has been carried out 
in a general covariant gauge. The results are provided as a polynomial of the clover parameter 
csw, and are tabulated for 10 popular sets of the Symanzik coefficients (Plaquette, Tree-level 
Symanzik, Iwasaki, TILW and DBW2 action). 

We also study the 0{a^) corrections to matrix elements of fermion bilinear operators that have 
the form ^F^, where Y denotes all possible distinct products of Dirac matrices. These correc- 
tion terms are essential ingredients for improving, to 0{a^), the matrix elements of the fermion 
operators. 

Our results are applicable also to the case of twisted mass fermions. 

A longer write-up of this work, including non-perturbative results, is in preparation together with 
V. Gimenez, V. Lubicz and D. Palao [0]. 
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1. Introduction 

Over the years, many efforts have been made for ^(a^) improvement in lattice observables, 
which in many cases is automatic by virtue of symmetries of the fermion action. According to 
Symanzik's program one can improve the action by adding irrelevant operators. Also, in the 
twisted mass formulation of QCD ^ at maximal twist, certain observables are ^(a^) improved, 
by symmetry considerations. 

So far, in the literature there appear two kinds of perturbative evaluations pertaining to the 
fermion propagator and bilinears of the form *Fr*P (F denotes all possible distinct products of 
Dirac matrices). On the one hand, there are 1-loop computations for ^{a^) corrections, with an 
arbitrary fermion mass [Q, |5|]. On the other hand, there are 2-loop calculations at level, for 

massless fermions 1-loop computations of G{cp-) corrections did not exist to date; indeed they 
present some novel difficulties. In particular, extending ^(a") calculations up to ^(a^) does not 
bring in any novel types of singularities. For instance, terms which were convergent to &(c^) may 
now develop an infrared (IR) logarithmic singularity at worst in 4 dimensions and the way to treat 
such singularities is well known. In most of the cases, e.g. for m = 0, terms which were already 
IR divergent to will not contribute to ^(a^), by parity of loop integration. On the contrary, 

the IR singularities encountered at 0'{a^) are present even in 6 dimensions, making their extraction 
more delicate. 



2. Description of the calculation 

Our calculation is performed for clover fermions, keeping the coefficient csw as a free pa- 
rameter. The action describing Nf flavors of degenerate clover (SW) fermions is given in Ref. |^]. 
We work with massless fermions (/Mo = 0), which simplifies the algebraic expressions, but at the 
same time requires special treatment for the IR singularities. By taking itIq = 0, our calculation and 
results are identical also for the twisted mass action in the chiral limit. 

For the gluon part we employ the Symanzik improved action, involving Wilson loops with 4 
and 6 links; for the Symanzik coefficients, a, multiplying each Wilson loop, we choose 10 sets of 
values that are widely used in numerical simulations; these are tabulated in Ref. [[I]]. 

The Feynman diagrams that enter this computation are shown in Fig. 1; diagrams 1 and 2 
contribute to the fermion propagator, while diagram 3 is relevant to the bilinears' improvement. 



1 2 3 

Fig. 1: Diagrams contributing to the improvement of the propagator (1,2) and the bilinears (3). 



A wavy (solid) line represents gluons (fermions). A cross denotes an insertion of F (Eq. (4-.1 )). 



For the algebraic operations involved in manipulating lattice Feynman diagrams, we make use 
of our symbolic package in Mathematica. Next, we briefly describe the required steps: 

• The evaluation of each diagram starts with the contraction among vertices, which is performed 
automatically once the vertices and the topology of the diagram are specified. The outcome of the 
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contraction is a preliminary expression for the diagram under study; there follow simplifications 
of the color dependence, Dirac matrices and tensor structures. We also fully exploit symmetries of 
the theory to limit the proliferation of the algebraic expressions. 

• The above simplifications are followed by the extraction of all functional dependencies on the ex- 
ternal momentum p (divergent, convergent terms) and the lattice spacing (terms of order a\ a^). 
The convergent terms can be treated by naive Taylor expansion in to the desired order. On the 
contrary, the isolation of the logarithms and non-Lorentz invariant terms is achieved as described 
below. As a first task we want to reduce the number of infrared divergent integrals to a minimal set. 
To do this, we use iteratively two kinds of subtractions for the propagator, so that all primitively 
divergent integrals (initially depending on the fermion and the Symanzik propagator) are expressed 
in terms of the Wilson propagator. The subtraction for the gluon propagator reads 

D(^) =Dp,,,(^)+Dp,,,(^)(Dp;i^(^)-D-i(^))D(^) (2.1) 

^j;(^) = ^-(l-^)|f, ^.=2sin(f), e = m (2.2) 

where D^^ is the 4x4 Symanzik propagator. The matrix {D~i^^{q) — (q)), which is ^{q^), is 
independent of the gauge parameter, A , and can be obtained in closed form, as a polynomial in . 

• The most laborious part is the computation of the divergent terms, which is performed in a nonin- 
teger number of dimensions D > 4. Ultraviolet divergences are explicitly isolated a la Zimmermann 
and evaluated as in the continuum. The remainders are D-dimensional, parameter-free, zero mo- 
mentum lattice integrals which can be recast in terms of Bessel functions, and finally expressed 
as sums of a pole part plus numerical constants. We analytically evaluated an extensive basis of 
superficially divergent loop integrals, which is presented in Ref. [|l|]. A few of these integrals are 
very demanding, because they must be evaluated to two further orders in a, beyond the order at 
which an IR divergence initially sets in. As a consequence, their evaluation requires going to D > 6 
dimensions, with due care to take into account all possible sources of i?{a^) corrections. These 
integrals form a sufficient basis for all integrals which can appear in any i?{a^) 1-loop calculation. 

• The required numerical integrations over loop momenta are performed by highly optimized For- 
tran programs; these are generated by our Mathematica 'integrator' routine. Each integral is ex- 
pressed as a sum over the discrete Brillouin zone of finite lattices, with varying size L (L^ < 128"^), 
and evaluated for all values of the Symanzik coefficients which we considered. 

• The last part of the evaluation is the extrapolation of the numerical results to infinite lattice 
size. This procedure entails a systematic error, which is reliably estimated, using a sophisticated 
inference technique; for 1-loop quantities we expect a fractional error smaller than 10^^. 

3. Correction to the fermion propagator 

The 1-loop corrections to the fermion propagator arise from the evaluation of diagrams 1 and 
2 in Fig. 1. Capitani et al. have calculated the first order terms in the lattice spacing for massive 
Wilson fermions and the Plaquette action for gluons. We carried out this calculation beyond the 
first order correction, taking into account all terms up to &{cp-) and considering a general Symanzik 
improved gluon action. Our results, to ^(a^), are in perfect agreement with those of Ref. |^. 
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The following equation is the total expression for the inverse propagator S^^ as a function of 
the external momentum p, the coupling constant g, the number of colors A^, the clover coefficient 
cswi and the gauge parameter A. The quantities e^''^^ appearing in our results for S^' are numerical 
coefficients depending on the Symanzik parameters, calculated for each action we have considered; 
they are tabulated in Ref. [[T]]. In Eq. (3^) and Eqs. (4^), we present the values of g^''^) for 
the Plaquette and Iwasaki actions (top and bottom numbers, respectively); only 5 decimal points 
are shown, due to lack of space 
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We define f = g^Cp /{I6n'^), Cp = {N^ - 1)/{2N), C2 = ci - C2 - C3 and / = j^p^ ; the spe- 
cific values A = 1 (A = 0) correspond to the Feynman (Landau) gauge. We observe that the tff{a^) 
logarithms are independent of the Symanzik coefficients c,- ; on the contrary iff{a^) logarithms have 
a mild dependence on c, . 
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4. Improved operators 

In the context of this work we also compute the contributions up to &{a^) to the forward matrix 
elements of local fermion operators that have the form *Fr*P. Y corresponds to the following set of 
products of the Dirac matrices 

r=l (scalar), 7^ (pseudoscalar), 7^ (vector), 7^7^ (axial), ^7^[7^,7v] (tensor) (4.1) 

The i?'{a^) correction terms are derived from the evaluation of diagram 3 shown in Fig. 1. One 
may improve the local bilinears by the addition of higher-dimension operators 

i=\ i=l 
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where the first term is the unimproved operator and are operators with the same 

symmetries as the original ones, but with dimension higher by one (two) units. To achieve &{a^) 
improvement, we must choose kj and kj appropriately, in order to cancel out all 0'{a^ ,a^) contri- 
butions in matrix elements. 

In the rest of this section we show our results for the 1-loop corrections to the amputated 2- 
point Green's function (diagram 3 of Fig. 1), at momentum p : A^{p) = (*P (*Fr*P) The 
values of all the Symanzik dependent coefficients, Es, Ep, Ey, £a, £t, with their systematic errors, 
can be found in Ref. We begin with the 6{c?-) corrected expressions for A'^(p) and A^{p)', 
including the tree-level term, we obtain 
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A^(7?) is free of i^(a^) terms and all contributions linear in csw vanish. The values of the numerical 
coefficients £5 and £p for the Wilson and Iwasaki gluon action are 
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The ^{a^) corrected expressions for (p), A^{p) and A^ {p) are very complicated, in the 
sense that there is a variety of momentum contributions and therefore many Symanzik dependent 

(4.7). In fact, we relegate our result for A^{p) to the 



coefficients, as can be seen from Eqs. 
longer write up We also list the coefficients £y and Sa for the Wilson and Iwasaki actions in 
Eq. (ph. 
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5. Conclusions 

Our results show clearly that i^{a^) effects are quite pronounced in the Green's functions we 
have considered. The &{a^) contributions which we have calculated can be directly used in order 
to construct improved operators, bringing the chiral limit within reach. Possible follow-ups to the 
present work include: 

• Extending to the case of nonzero renormalized mass. 

• Improvement of higher-dimension bilinear operators, such as those involved in hadronic form 
factors, and of 4-fermi operators. 

A comparison with non-perturbative estimates of matrix elements, coming from numerical simula- 
tions, will be presented in Ref. [[T]]. 
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